The Dynamics of ID Quantum Spin Systems Can Be Approximated Efficiently 
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In this Letter we show that an arbitrarily good approximation to the propagator e**^ for a ID lattice of 
n quantum spins with hamiltonian H may be obtained with polynomial computational resources in n and 
the error e, and exponential resources in Our proof makes use of the finitely correlated state/matrix 
product state formalism exploited by numerical renormalisation group algorithms like the density matrix 
renormalisation group. There are two immediate consequences of this result. The first is that Vidal's time- 
dependent density matrix renormalisation group will require only polynomial resources to simulate ID 
quantum spin systems for logarithmic \t\. The second consequence is that continuous-time ID quantum 
circuits with logarithmic \t\ can be simulated efficiently on a classical computer, despite the fact that, after 
discretisation, such circuits are of polynomial depth. 

PACS numbers: 75.10.Pq, 03.67.-a, 75.40.Mg 



The kinematics and dynamics of quantum lattice sys- 
tems are strongly constrained by a key physical require- 
ment, namely, the locality of interactions. For example, 
consider a collection of n distinguishable spin-i systems 
which interact according to nearest- neighbour interactions: 
a counting argument quickly reveals that such local hamil- 
tonians occupy an extremely small fraction of the space of 
general hamiltonians. Thus, intuitively, we would expect 
this nongeneric constraint would manifest itself strongly in 
the structure of the eigenvalues and eigenstates for local 
hamiltonians. This is indeed the case, but it is still far from 
obvious exactly how to best quantify this constraint. 

A number of methods to systematically quantify the 
eigenstates and matrix functions of local hamiltonians have 
been developed. Perhaps the most successful scheme in re- 
cent years has been the technology of finitely correlated 
quantum states (FCS) |1, 2]. (Finitely correlated states are 
also known as matrix product states (MPS) in one dimen- 
sion and tensor product states or pro jected entangled-pair 
states in two and higher dimensions |3]. The key feature of 
a finitely correlated state is that, as the name suggests, sep- 
arated regions are weakly correlated. In addition, any state 
which does not exhibit too much correlation between sep- 
arated subsystems can be well approximated by a finitely 
correlated state |4], |22]. 

Finitely correlated states are nothing more than a con- 
venient representation for vectors in tensor-product hilbert 
spaces. However, the utility of this particular representa- 
tion is that for those states with bounded or limited corre- 
lations it is often extremely efficient (in n) to extract local 
properties, such as expectation values of local operators. 

The utility of the FCS representation as a means to cal- 
culate local properties of ID quantum lattice systems has 
been spectacularly demonstrated by the development of the 
density matrix renormalisation group (DMRG). (See |5] 
and references therein for a description of the DMRG and 
related algorithms.) The DMRG provides an apparently ef- 
ficient computational recipe to obtain an approximation to 
the ground state and low-energy eigenstates for ID quan- 



tum lattice systems as FCS vectors. The DMRG is an ex- 
tremely flexible method and has been recently extended to 
apply to a diverse number of situations, such as the calcu- 
lation of short-time dynamics fl0l, dissipation |8, 9], dis- 
ordered systems 1 10], eigenstates with definite momentum 
lUI, and, recently, higher dimensions [3]. 

Perhaps one of the most exciting recent results in the 
study of the DMRG has been the development of an al- 
gorithm to simulate the real-time dynamics of ID quantum 
spin systems The efficiency of this algorithm is 

predicated on the condition that the dynamics of the spin 
system do not create too much long-range quantum en- 
tanglement 10, 0|- While it appears that, in practice, this 
condition is always fulfilled for small times, it is currently 
unclear if it is true for all ID local quantum spin systems. 

There are at least two reasons why it is interesting to 
study the theoretical worst-case computational costs of the 
time-dependent DMRG. The first is that the computational 
complexity of the DMRG and related algorithms is cur- 
rently unknown except for when applied to a handful of 
singular integrable models 1 1 3 ] . An assessment of the theo- 
retical worst-case computational complexity of the DMRG 
in any other circumstance would allow one to certify a 
priori the accuracy of the DMRG versus computational 
cost. The second reason is related to computational power 
of quantum computers (see |14] for a detailed descrip- 
tion of quantum computation and a number of quantum al- 
gorithms including quantum simulation algorithms). The 
time-dependent DMRG provides a way to simulate quan- 
tum computers running quantum algorithms. A careful the- 
oretical worst-case complexity analysis would potentially 
give us an insight into what quantum computations can and 
can't be simulated efficiently on a classical computer. 

The most naive way to study the computational com- 
plexity of the time-dependent DMRG is to directly study 
the storage costs of representing the propagator e'*^ when 
it is approximated by the Lie-Trotter expansion e'*^ ~ 



(e*'" e*™ )™, for some large m 12311 . Unfortunately, 
however, a careful analysis of the error scaling with m for 
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FIG. 1 : Decomposing the propagator: for a given time t the prop- 
agator e'*^ may be written as a product of e^'^^-^+^s) and V{t), 
where Hy^ (respectively, Hg) is the hamiltonian with interaction 
terms only between spins in region A (respectively, between spins 
in region B) and V{t) is a unitary operator which "patches up" 
the error due to approximating e**^ with e**(^^+^B). Because 
of the UV cutoff given by the lattice, information propagation in 
quantum spin systems is limited by an effective "speed of light". 
Hence the unitary V{t) interacts spins in A with spins in B suc- 
cessively weaker further from the boundary between A and B. 
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FIG. 2: The approximate quantum cellular automata decomposi- 
tion for e'*^ for time t. Each of the unitary operators U and V 
acts on at most co\t\ spins, where cg is a constant. 



this representation shows that the worst-case storage cost 
might be exponential, even for |t| which scales as a con- 
stant with n. This, in turn, implies that quantum circuits 
simulating ID dynamics using the Lie-Trotter expansion 
have a depth that scales at least linearly with n. If we 
want to simulate such methods efficiently we need a more 
sophisticated technique to obtain a representation for the 
propagator. 

In this Letter we show that the propagator e'*^ for an 
arbitrary ID quantum system is well- approximated by a 
finitely correlated state vector in the Hilbert space of opera- 
tors using resources that scale polynomially with n and ex- 
ponentially with Because our proof is constructive we 
obtain an efficient algorithm, closely related to the DMRG, 
to obtain this representation. 

The argument we describe in this Letter can be under- 
stood by appealing to the following physical intuition. The 
dynamics of any ID quantum spin system are constrained 
by the ultraviolet cutoff induced by the lattice spacing be- 
tween the spins. This cutoff induces a bound on the speed 
at which information can propagate in such systems, an in- 
tuition which is precisely quantified by the Lieb-Robinson 
bound 1 15]. We exploit this bound on information propaga- 
tion to provide two different decompositions of the propa- 
gator e'*^. In the first decomposition we partition the chain 
into two contiguous pieces A and B and approximate the 
dynamics e'*^ by e**(^-*+^^\ where Ha is the hamilto- 
nian which includes interaction terms only between spins 



in A (and similarly for Hb)- Obviously this approximation 
is not exact: at the cut point between A and B there will 
be substantial errors. However, as a consequence of the 
bound on information propagation, the difference between 
the way e'*^ and ^-^KHa+Hb) ^^.j spins far from the 
cut will become small (information about the cut doesn't 
have time to propagate too far away). We can patch up the 
difference between the two propagators by introducing a 
new unitary operator V{t) which acts across the boundary: 
^itH ^ ^it{HA+HB)Y(^^y The Lieb-Robinson bound then 
tells us that V{t) acts progressively weaker on spins far 
from the cut (see Fig. ^ for discussion and illustration of 
this). This, in turn, allows us to approximate V{t) with a 
unitary V'{t) which acts nontrivially only on a finite num- 
ber of spins i7 around the boundary. If |0| is bigger than 
\t\ then this approximation improves exponentially fast in 

We obtain the second decomposition for e'*^ by first fix- 
ing t and then moving along the chain C and introducing a 
cut after approximately |t| spins (which we call Ai) and 
then patching it up with a unitary Vi (t) which acts nontriv- 
ially only on a finite number of spins across the boundary 
(this is the unitary constructed in the previous paragraph). 
Thus we have e**^ ~ yi(t)e'*^^i O e^*^^\^i. Then we 
recursively apply this procedure to e**^'^\'^i until we arrive 
at the approximate decomposition illustrated in Fig. [J] 

We will, for the sake of clarity, describe our results 
mainly for a finite chain C of n distinguishable spin- 1/2 
particles. The family H of local hamiltonians we focus 
on (which implicitly depends on n) is defined hy H = 
Yl^=o ^3 ' where hj acts nontrivially only on spins j and 
j + We set the energy scale by assuming that \\hj\\ scales 
as a constant with n for all j = 0, 1, . . . , n — 1, where || • || 
denotes the operator norm. The interaction terms hj may 
depend on time: hj = hj{t). We can easily accommodate 
next-nearest neighbour interactions etc. by blocking sites 
and thinking of the blocks as new (larger) spins. However 
this can only be done a constant number of times: the qual- 
ity of our approximation will decrease exponentially with 
the number of such blockings. We do not assume transla- 
tional invariance. 

The crucial idea underlying our approach is that a good 
approximation to the propagator e**^ for a local ID quan- 
tum spin lattice system can be obtained and stored effi- 
ciently (i.e. with polynomial resources in n) with a clas- 
sical computer for |i| < clog(n), where c is some con- 
stant. The way we do this is to use a specific representation 
for the approximation, namely as a finitely correlated state 
vector. What we mean by this is that we represent an oper- 
ator W in the following fashion 

W = ^ A"" • • • A""- V"" 0-"^ ® • • • (g) cj""-^ , 

(1) 

where Q„ = {0,1,2,3}''", A"" (respectively, A""-i) 
are a collection of four row vectors of size (respectively. 



3 



four column vectors of size C„_i), A"^ are four Cj x Dj 
sized matrices, for 1 < j < — 1, and cr" is the vector 
of Pauli operators. Note that C^+i = Dj. The dimensions 
Cj and Dj are called the auxiliary dimensions for site j. 
It is clear that if the sizes of the auxiliary dimensions are 
bounded by polynomials in n, i.e. if Cj < poly(n) and 
Dj < poly(n), then the operator W can be stored with 
polynomial resources in n. Also note that all operators can 
be represented exactly as in Eq. ([B by taking the auxiliary 
dimensions to be large enough: Cj = Dj = 2" suffices 

We begin by showing how to obtain an approximate de- 
composition of the propagator e**^ as a product 



^n/|n|-l 



itH 




Vn'(t)]+Oie), (2) 



where we have two partitions Vi and 7^2 of the chain C 
into n/jr^l (respectively + 1) contiguous blocks of 

< \Q\ spins. The first set is denoted Vi = {^j}- The 
second set 1^2 = {^'j} is the set of blocks which are just 
translates of those in Vi by \^\/2 sites fl^. The oper- 
ators Un- (respectively, Vq'.) are unitary operators which 
act nontrivially only on Qj (respectively, Q'j). We call 
such a decomposition an approximate quantum cellular 
automata decomposition (or, simply, an e-QCA decom- 
position) because e**^ is exactly a Margolus-partitioned 
QCA update rule (see fvf\ for a description of QCAs). 
Then we show that an e-QCA decomposition implies that 
^itH well-approximated by a PCS vector with auxil- 
iary dimension 21^'. To reduce the error to e we require 
I^^I ^ Co[t[ + Ci log(n/e), where Cq and Ci are constants. 
This decomposition is illustrated in Fig. [J] 
Consider the unitary operator 



V{t) = (e~'*^^\^ 



(3) 



As we described in the introduction, for small and for 
sites far enough away from the boundary dh between C \ 
A and A, this operator ought to be close to the identity. 
Therefore we argue that, as an operator, V{t) ought to be 
expressible as V{t) ~ Ic\o ^(^)- 

To quantify this statement we study the differential equa- 
tion that V{t) satisfies: 



dV{t) 
dt 



JtH 



le 



iVityl'ihj), (4) 



where hj is the interaction term that bridges the left- and 



its 



right-hand side of the chain, and Tf{A) = j4e 
Thus we see that V{t) is generated by time-dependent uni- 
tary dynamics due to the effective hamiltonian L{t) = 
T^[hi), and we write 

=Te''/o^(^)''«, (5) 

where T denotes time-ordering. 

The time-dependent effective hamiltonian L{t), and 
hence V{t), acts nontrivially on all of the sites in the 



chain. We now show that, for small L{t) may be 
well-approximated by an operator which acts nontrivially 
on only a handful of sites near the boundary dh between 
the left- and right-hand sides of the chain. To do this, we 
construct the following approximation to L{t): 



L'{t) 



{hi 



(6) 



where i7 is a contiguous block of sites centred on dA, and 
contains only those interaction terms hj in H which 
interact only spins contained in 17. We now show that, 
for small enough |t| and large enough 17 containing Hj, 
\\L{t) — L'{t)\\ < e, for some prespecified e. 

To show that L'(t) is a good approximation to L(t) we 
must establish that \\t^ (hj) —t^''' {hj)\\ is small. Abound 
on such a quan tity is known as a Lieb-Robinson bound 
flS*, Tl in ll^] (see in for a simple direct proof). The 
strongest (and easiest to prove) such bound reads 



tj.\n\ 



i=\n\ 

where 5, = ||/ij||2'||/i||', \\h\\ = 
ui are constants. Thus we find 

\\L{t)-L' 



max,- \\hj\\, and k, and 



(7) 



In this way we see that we can reduce the operator norm 
difference between L[t) and L'{t) exponentially fast in the 
size 1 17 1 of the region 17. 

We now define a new unitary operator V'{t) — which 
is meant to approximate V{t) — as the unitary operator 
generated by the time-dependent hamiltonian L'(t): 



dV'jt) 
dt 



iV'{t)L'{t). 



(8) 



Because L'{t) acts nontrivially only on 17, V'{t) is a uni- 
tary operator which acts nontrivially only on 17 and it acts 
as an identity elsewhere. In order to see how accurate V{t) 
is as an approximation to V'{t) we now bound the error 
\\V(t)-V'(t)\\. 

To show that V{t) and V'{t) are close for some time 
period we integrate the differential equations © and (|8j. 
We do this by making use of the Lie-Trotter expansion 



V{t) = lim 



J=0 



V'{t) = lim TT e'^' 



(9) 



(10) 



3=0 



applying the triangle inequality several times, and taking 
the limit m — > cxo. This gives us the fundamental estimate 

\\V{t)-V'it)\\ < r\\Lis)-L'{s)\\ds. (11) 



Substituting (Ql and redefining constants gives us 

\\V{t) - V'{t)\\ < we'l'le-^l'^l, (12) 
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where uj, k, and /i are constants independent of n. 

Our final result is now the following. Rearranging (|3j 
and using the estimate (I12t gives us 

e''" = (e^*^'^\^ (8) e'*^^)y'(t) + e, (13) 

where V'{t) acts nontrivially only on a contiguous block 
of spins of size Iterating this procedure by cutting A 
into two pieces etc. give us the final e-QCA decomposition 
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(n/\n\-i \ /n/\n\ \ 

® \ ®J^^'Mj+Oie), (14) 

where = O(co|t| + Cilog(n/e)), for some constants 
Co and Ci . 

It is now relatively straightforward to show that an e- 
QCA decomposition gives rise to an efficient ECS vector 
representation once we recognise that the expression 

/n/\n\-i \ 

U{t) =y® Un,{t)j (15) 

is a ECS vector with auxiliary dimension 2'^^' The way to 
see this is to note that, in the standard operator basis, 

n/|n|-l / \ 

n E 

3=0 \k,eQ|n| / (16) 

k 

which is a ECS vector representation with maximum aux- 
iliary dimension 2'^'. 
Given that 

/n/\n\ \ 
V{t)= \ ®Vn'^{t)j (17) 

is also expressible exactly as a ECS with maximum aux- 
iliary dimension 21^' and using the result that the product 
of two ECS vector operators with maximum auxiliary di- 
mensions Di and D2 admits a ECS vector expression with 
maximum auxiliary dimensions D1D2 we obtain the final 
result that our approximation V'{t) = lA{t)V{t) to V{t) 
is expressible exactly as a ECS with maximum auxiliary 
dimension < 2^'^'. If |t| scales as a constant, or logarith- 
mically, with n then we learn that the ECS representation 
requires only polynomial resources in n. 

We have shown how the propagator e**^ for a ID sys- 
tem of quantum spins can be efficiently obtained and rep- 
resented using a classical computer. There are many con- 
sequences of this representation. The first and most ob- 
vious is that the time-dependent DMRG will, in the theo- 
retical worst case, require only polynomial computational 
resources to simulate time evolution for constant times. 
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